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The role of exterior Weyl fluids on compact 
stellar structures in Randall-Sundrum gravity 
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In the context of the Randall-Sundrum braneworld, the minimal geometric deformation approach 
(MGD) is used to generate a new physically acceptable interior solution to Einstein's field equations 
for a spherically symmetric compact distribution. This new solution is used to elucidate the role 
of exterior Weyl stresses from bulk gravitons on compact stellar distributions. We found strong 
evidences showing that the exterior dark radiation IA + always increases both the pressure and the 
compactness of stellar structures, and that the exterior "dark pressure" V + always reduces them. 
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I. INTRODUCTION 

General Relativity (GR) represents, without any 
doubt, one of the most important achievement of Physics. 
The predictions made for this theory, like the perihelion 
shift of Mercury, light deflexion, gravitational red shift, 
gravitational lens, time delay, etc, have given to it the 
honor which deserves as one of the fundamental theories 
of Physics (for an excelent review on experimental tests of 
GR see [lj and references therein). On the other hand, 
due to the great technological advance during the past 
years, today we have increasingly powerful instruments 
which allow to obtain accurate measures associate to the 
evoution of compact objects, which represent excellent 
laboratories for the study of gravity in the strong regime 
(see for instance 0). This not only serves to test the 
theory like never before, but also leaves GR as the only 
reliable gravitational theory to be used in the analysis of 
phenomena occurred in the strong field regime. Likewise, 
the ability of observing increasingly distant objects deep 
in the universe, and thus with a great gravitational red 
shift, leads inevitably to the conclusion that only using 
GR we can obtain an adequate analysis of these phenom- 
ena [3], [J]. Furthemore, with the recent results shown 
by PLANCK [fj, which improve greatly the previous by 
WMAP 0, we can assure that the cosmological mod- 
els based in GR enjoys a well-deserved and well-earned 
prestige. 

Despite the above facts, there are some fundamen- 
tal questions associated to the gravitational interaction 
which GR cannot answer satisfatory. This can be broadly 
grouped in two fundamental issues, which most likely are 
closely related: 1) The inability of GR to explain satis- 
factorily the dark matter [7( and dark energy problem 
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without the need of introducing some kind of unknown 
matter-energy to reconcile what predicts GR with the 
observed, namely, galactic rotation curves and acceler- 
ated expansion of the universe. 2) The impossibility, so 
far, to reconcile GR with the Standard Model of parti- 
cle physics, or equivalcntly, the inability to quantize GR. 
This has strongly motivated the searching of a gravita- 
tional theory beyond GR that helps to explain satisfac- 
torily part of the problems described above. If the new 
theory is a consistent quantum theory, this should lead 
to a generalization of GR at low energy, being likely this 
extension of GR at low energy which could accounts for 
the dark matter and dark energy problems. If the new 
theory is not a consistent quantum theory for gravity, this 
should also contain GR in a suitable limit, and somehow 
show greater tolerance to its quantum description. 

Extra-dimensional theories, which arc mostly inspired 
by String/M-theory, are among the theories that lead to 
modifications to GR. One of these extra-dimensional the- 
ories is the Braneworld (BW) proposed by Randall and 
Sundrum (RS) Q which has been largely studied and 
which explains, so pretty straightforward, one of the fun- 
damental problems of Physics, i.e. the hierarchy prob- 
lem (see also the ADD model Q). This theory reduces 
the fundamental scale to the weak scale by considering 
extra-dimensional effects, thus explaining the weakness 
of gravity relative to the other forces. Because of this, its 
study and impact on GR is fully justified and is of great 
importance [10( , jllj , although not yet found experimen- 
tal evidence to support it [12(, (l3j |. 

Despite the great efforts made in recent years and the 
great understanding of the theory of RS braneworld, we 
are still far from fully understanding its impact on grav- 
ity, mainly in self gravitational systems. Although we 
have a covariant approach that is useful to study many 
fundamental aspects of the theory [lj] , there are certain 
key issues that remain unresolved. One of these problems 
is the necessity to clarify the existence of black holes so- 
lutions in RS. In order this theory can be considered a 
physically consistent theory, it must support, at least, 



black holes solutions. Today we have some evidences 
indicating the existence of black holes solutions in RS 
theory [15l - ll7| . but an exact solution remains unknown 
so far. Indeed, solving the full five-dimensional Einstein 
field equations, something which would be tremendously 
useful, has proven to be an extremely complicated prob- 
lem (see, e.g. [la 13; an d references therein). Another 
aspect worth noting is that, even knowing an exact solu- 
tion, it would produce different solutions in our observed 
universe, since there are many ways to embed a four- 
dimensional brane in the five-dimensional bulk. There 
are indeed many different ways of embedding, e.g. iso- 
metric, conformal, imbeddings, rigid, global, local, an- 
alytic etc. [201. While it is true that Z2 symmetry is 
popular when considering this point, it is far from being 
a problem widely studied. 

Since the complete solution (bulk plus brane) remains 
unknown so far, finding exact solutions, or at least physi- 
cally acceptable solutions to effective Einstein field equa- 
tions in the brane, might be certainly a good guide to 
clarify some aspects of the five-dimensional geometry and 
the way our observed universe is embedded in it. Once 
a solution is found, we could use Campbcll-Magaard 
theorems [2l], (22| to extend the brane solution through 
the bulk, locally at least (To see some consequences of 
Campbell-Magaard theorems to general relativity, see 
Ref [2J,[2J]). To accomplish this, first of all we need to 
start with a rigorous study of the effective Einstein field 
equations in the four-dimensional brane. A study that, 
among other things, clarify the role of five-dimensional 
effects on the effective four-dimensional field equations, 
and above all, a study to consolidate a general method- 
ology based on a critical basic requirement: regain GR 
at low energies. Fortunately, the RS theory has a free 
parameter, the tension of the brane, which can be used 
to control this aspect of extreme importance, and which 
is a nontrivial problem [251 ] . In fact, this fundamental 
requirement is the basis of the minimal geometric defor- 
mation approach (MGD) |26[, which has allowed, among 
other things, to generate physically acceptable interior 
solutions for stellar systems [27|, physically acceptable 
exact interior solutions (28[, to solve the tidally charged 
exterior solution found in Ref [29| in terms of the ADM 
mass and to study (micro) black hole solutions [30, bLL| . 

In the RS BW theory, there are two fields filling our 
four-dimensional vaccum, namely, the dark radiation U 
and dark pressure V , which have an extradimensional ori- 
gin, and whose effects on stellar structures is not well un- 
derstood so far. In this paper the MGD approach is used 
to generate a new physically acceptable interior solution 
to four-dimensional effective Einstein's field equations for 
a spherically symmetric compact distribution, which is 
used to elucidate the role of exterior Weyl stresses from 
bulk gravitons on compact stellar distribution. This pa- 
per is organized as follows. In Section II the Einstein 
field equations in the brane for a spherically symmetric 
and static distribution of density p and pressure p is re- 
minded, as well as the MGD approach. In Section III a 



new physically acceptable stellar interior solution to Ein- 
stein's fields equations is generated by breaking a well 
known general relativistic perfect fluid solution through 
the MGD approach. In Section IV the consequences of 
dark pressure V and dark radiation U on stellar struc- 
ture will be analized through the matching conditions by 
using two different exterior solutions. In the last section 
the conclusions are presented. 



II. FIELD EQUATIONS AND THE MINIMAL 
GEOMETRIC DEFORMATION APPROACH. 

In the context of the braneworld, the five-dimensional 
gravity produces a modification on Einstein's field equa- 
tions in our (3+1) observed universe, the so-called brane. 
These modifications can be seen through the energy- 
momentum tensor, which now has new terms carrying 
five-dimensional consequences onto the brane: 
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where a is the brane tension, with S^ and £ M „ the high- 
energy and non-local (from the point of view of a brane 
observer) corrections respectively, and J- ' ^ a term which 
depends on all stresses in the bulk but the cosmological 
constant. In this paper, only the cosmological constant 
will be considered in the bulk, hence T^v = 0, which 
implies there will be no exchange of energy between the 
bulk and the brane, and therefore V T M „ = 0. 

Using the line element in Schwarzschild-like coordi- 
nates 
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the effective equations in the brane for a spherically sym- 
metric and static distribution, with Weyl stresses in the 
interior, can be written as 12611 
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where /' = df/dr and k 2 = 8ir. The general relativity 
is regained when ct _1 — > and Eq. ([5]) becomes a linear 
combination of Eqs. (J3j)- dSJ) . 

Despite the fact that Eqs.©-© represent an indefi- 
nite system of equations in the brane, an open problem 
for which the solution requires more information of the 
bulk geometry and a better understanding of how our 4D 
spacetime is embedded in the bulk, it is possible to gen- 
erate the braneworld version of every general rclativistic 
solution through the MGD approach (26[. In this ap- 
proach, from the point of view of a brane observer, the 
five-dimensional gravity produces a geometric deforma- 
tion in the radial metric component given by 

" pdr + Geometric Deformation . 

" v ' 

(9) 

When a solution of 4D Einstein's field equations is con- 
sidered as a possible solution to the braneworld system, 
given by Eqs. ©-(JB]), the geometric deformation in Eq. 
((9) is minimal, and it is given by 
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As it is shown by Eq. (|T0) . the geometric deformation 
f*(r) satisfies f*(r) Js 0, hence it always reduces the 
effective interior mass, as it is seen further below in Eqs. 
(|T8)l and ([19]). On the other hand, Eq. ((TO) represents 
a minimal geometric deformation in the sense that all 
sources of the geometric deformation have been removed 
except those produced by the density and pressure, which 
are always present in a stellar distribution. However, 
there is an even minimal deformation in the case of a dust 
cloud p = 0, but this will not be considered in the present 
work. This geometric deformation /* (r) is the source of 
the anisotropy induced in the brane, whose explicit form 
may be found through Eq. (J4), leading to 

K a \ a \rr 



(12) 

It is clear that this minimal deformation will produce a 
minimal anisotropy onto the brane. 



III. A STELLAR SOLUTION. 

The Hcintzmann solution [32[ is a well known spher- 
ically symmetric stellar solution for a perfect fluid in 



general relativity. This solution satisfies all the elemen- 
tary criteria for physical acceptability that a stellar solu- 
tion must satisfy, namely, regular at the origin, pressure 
and density defined positive, well defined mass and ra- 
dius, monotonic decrease of the density and pressure with 
increasing radius, dominant energy condition satisfied, 
subliminal sound speed, etc. Because of the above, the 
Hcintzmann solution represents a very attractive candi- 
date to be considered in the braneworld context through 
the minimal geometric deformation approach. Hence 
some interesting aspect of the five-dimensional gravity 
consequences on stellar structures might be elucidated. 

Let us start by considering the Hcintzmann solution 
for a perfect fluid in general relativity 
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where A, B and C are constants to be determined by 
matching conditions. In general relativity, all these con- 
stants have specific values. Indeed, they may be written 
in terms of the compactness of the distribution, that is, 
in terms of M/R, with M and R the mass and radius 
of the distribution, which are free parameters satisfying 
the constraint M/R < 4/9. However, as it is well known, 
in the braneworld scenario the matching conditions are 
modified, consequently there are five-dimensional effects 
on these constants which must be considered. 

In general relativity the second fundamental form, 
which leads to p(r) \ r =R = at the stellar surface r = R, 
produces 
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We will keep the physical pressure vanishing on the sur- 
face, even though this condition may be dropped in the 
braneworld scenario [33[ , [3J] . 

From the point of view of a brane observer, the geo- 
metric deformation f*(r) produced by five-dimensional 
effects "breaks" the perfect fluid solution represented by 
Eqs. (|13p -(fT6 ) . introducing thus imperfect fluid effects 
through the braneworld solution for the geometric func- 
tion A(r), which is obtained using Eqs. (TT3) . (|15[) and 
dHU in Eq. ©, leading to 
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where the interior mass function rh is given by 



m(r) = m(r) - - f*(r), 



(19) 



with f*(r) the minimal geometric deformation for the 
Heintzmann solution, given by Eq. (|10[) . whose explicit 
form is obtained using Eqs. (fT3]) . ([T5| and (fT6|) in Eq. 
(fTU) . hence 
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The function m(r) in Eq. (|19[) is the general relativity 
interior mass function, given by the standard form 
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hence the total general relativity mass is obtained 
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On the other hand, it can be shown by Eqs. (jTTj) and 
(J2"D"|) that the geometric deformation f*(r) depends only 
on the parameter C, which has a well defined expression 
in terms of the compactness in general relativity. The 
general relativity expressions for A and C, which from 
now on will be called Aq and Cq respectively, are found 
in terms of M and R by the continuity of the metric at 
the stellar surface, where, of course, the Schwarzschild's 
exterior solution must be used. Thus, by considering the 
temporal component of the metric, we have 
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Hence using Eq. <(T7]) in Eq. ([22} and Eq. ([23) two simple 
expression relating Aq and Co with the compactness are 
found as shown below 
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the geometric deformation. A numerical analysis carried 
out inside the stellar distribution, for different densities, 
shows that the anisotropic stress V is proportional to 
the density: the most compact distribution undergoes a 
higher anisotropic effect. This behavior, which has al- 
ready been observed in the study of uniform stellar dis- 
tributions [23], can be easily explained in terms of the 
source of the anisotropy, which is the minimal geometric 
deformation undergone by the radial metric component, 
explicitly shown through Eq. ([TU)) . As the Heintzmann 
solution is a solution to AD Einstein's field equations, it 
removes all the non-local sources from the geometric de- 
formation f(y, p,p) in the generic expression given by Eq. 
([9"|h leaving only the high energy terms shown explicitly 
in Eq. (j 10[) , which are quadratic terms in the density and 
pressure. Hence the higher the density, the more geomet- 
ric deformation will be produced, and as a consequence 
the anisotropy induced will be higher for more compact 
distributions. On the other hand, when the scalar Weyl 
function U is analyzed, it is found that this function is 
more negative for more compact stellar objects. The fact 
that U is always negative means that high energy terms 
always dominate anisotropic terms, which are the two 
sources for U, as can be seen through Eq. ([5]). Next, us- 
ing two different exterior solutions having a Weyl fluid, 
the consequences of dark pressure V and dark radiation U 
on stellar structure will be analizcd through the matching 
conditions. 
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showing thus that the geometric deformation f*(r) in 
Eq. (|20p can be written in terms of the compactness of 
the distribution, therefore it may be used as a free pa- 
rameter to model different compact stellar distributions. 
Indeed, it is found that more compact distributions un- 
dergo a higher deformation due to five-dimensional ef- 
fects. In order to obtain the interior Weyl functions V 
and U, Eqs. (Q3J), ([15) and $jo$ along with Eq. (JTS) and 
Eq. (Q33) are used in Eq. (@| and Eq. ([5), leading to ex- 
pressions too large to be shown here, but proportional to 



IV. THE ROLE OF EXTERIOR WEYL FLUIDS 

It is well known that in general relativity the unique 
exterior solution for a spherically symmetric distribu- 
tion is the Schwarzschild metric. This situation changes 
dramatically in the five-dimensional braneworld sce- 
nario, where high energy corrections and the presence 
of the Weyl stresses imply that the exterior solution 
for a spherically symmetric distribution is no longer the 
Schwarzschild metric. Among all possible exterior so- 
lutions known to date, the one obtained by Dadhich, 
Maartens, Papadopoulos and Rezania (DMPR) in Ref. 
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represents the simplest generalization of the 
Schwarzschild exterior solution in the braneworld. 
Using our interior solution given by Eqs. (|13|) and (|18p 
along with the Reissner-Nordstrom-like solution Eq. 
([26]) in the matching conditions at the stellar surface 
r = R, we have 
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The tidal charge q in Eq. Q28) and Eq. <J29J) is ob- 
tained from the second fundamental form, which in the 
braneworld is given by 



£-£-™+£- 



PA+- 



1 



P% 



PrPr 



2 



4_V 
k 4 






but in this approach is reduced to 
f*(C) (v'(R) 1 \ 



PR 



Sn 



R 



R? 



k 4 



R 

a 



2_Ul 
k 4 a 



±K 

k 4 a ' 



thus using the DMPR solution given by Eq. 
(13TT) we found 
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which for our interior solution given by Eq. (|13[) leads to 
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The expression f*(C) = f*(R) is the minimal geometric 
deformation at the stellar surface r = R. Since the ge- 
ometric deformation f*(C) is always positive, it can be 
seen from Eq. ([33)) that the tidal Weyl charge q is always 
negative. The constants A4 and q are given in terms of C 
through Eqs. (|29p and ([33)) respectively, and C may be 
determined by Eq. ([28| if /I is kept as a free parameter, 
which can be used to find a physically acceptable model. 
On the other hand, the general relativistic value of C, 
named Co, can be seen by Eq. (|28p evalated at a^ 1 = 0, 
which leads to Eq. ([23)) . Hence comparing Eq. ([28)) with 
Eq. ([2"3"| , it is clear that the general relativistic values of 
C has been modified by bulk effects as 



C(a)=C + S(a). 
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In order to determine 8(a), Eq. ([34)) is used in Eq. (J28I) . 
hence 

A[l + (C a + 8)R^ = l- 2 -^ + ^L. (35) 

Now using Eqs. ([29)) and (|33)) in Eq. ([351 . and keeping 
linear terms in <r _1 , we found 
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As expected, the modification of C is proportional to the 
geometric deformation at the surface. 

At this stage, we have all the necessary tools needed 
to examine the braneworld consequences on the physical 
variables. For instance, to see five-dimensional conse- 
quences on the pressure p for different distributions, all 
we have to do is to use Eq. ([3"6")) in Eq. ([TB)) with differ- 
ent values of the positive parameter A (A < 1). As the 
compactness of the distribution is increased when A de- 
creases, we may obtain the five-dimensional consequences 
for different compact distributions. Figure 1 shows the 
behavior of the pressure in both the general relativity 
and braneworld case. In the braneworld case, two dif- 
ferent exterior solutions are considered. It can be seen 
that the pressure is increased by five-dimensional effects 
when the DMPR solution is used. The same result was 



found in Ref [27] for a uniform stellar distribution de- 
scribed by a braneworld version of the Schwarzschild's 
metric, where high energy corrections were considered 
along with Weyl stresses from bulk gravitons, which rep- 
resents a different result from other braneworld solution 
[35|, where only high energy modifications were consid- 
ered. The latter strongly suggests that, at least for the 
Schwarzschild's braneworld model, unlike the high energy 
effects, the Weyl stresses increases the pressure, and that 
its effects dominate over the high energy effects. How- 
ever, the specific role played by each Weyl function can- 
not be elucidated yet, unless one of them is turned off 
in order to see its consequences, but accomplishing this 
without breaking the interior braneworld solution is very 
complicated. The only way to see the individual effects 
of Weyl functions on stellar structure is through an ex- 
terior solution with one Weyl function off. Fortunately, 
there is a vacuum solution that complies with this fea- 
ture, namely, the vacuum braneworld solution foun d by 
Casadio, Fabbri and Mazzacurati (CFM) in Ref. [3£ 
given by 
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where r\ is a constant which measures deviation from 
the Schwarzschild's solution. Using our interior solution 
given by Eq. flTSJ) and Eq. (ISJ) with the CFM solution 
Eqs. ([3"7| - ([3"8"[) in the matching conditions at the stellar 
surface r = R, we have 



with the second fundamental form Eq. (|3Tj) leading to 
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The constants M. and 77 are obtained in terms of C by 
Eq. (g5]) and Eq. ([15)1. while C is determined through 
Eq. ([4"Tj) . where A is kept as a free parameter which 
again will be useful to modeling distributions with differ- 
ent compactness. Using Eq. (|42)) in Eq. (|4Tj) we obtain 
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Keeping linear terms in cr 1 and 77, Eq. (|44[) may be 
(42) written as 

I 
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Comparing this last expression with its general relativistic counterpart given in Eq. (|23p , it is clear that the general 
relativistic values of C has been modified by bulk effects by 



C(a)=C a + 5(a,i 1 ). 
Using Eq. (|40[) in Eq. (|45p and keeping linear terms in <5, we obtain 
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In order to find 77 = 77(a), Eq. (|42[) is used in the second 
fundamental form Eq. (|43[) . hence 
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From the last expression and since the geometric defor- 
mation /* > 0, it can be seen that 77 < 0, hence the 
exterior solution has a singular horizon at r = j^- [36( | . 
Using Eq. (|48[) in Eq. (|47|l and keeping linear terms in 
<7 -1 , after some algebraic manipulation and considering 
the zero order of Eq. (|4"5)l . which is given by Eq. ([2"3"]l . it 
is found that under the CFM exterior solution the devi- 
ation of C is given by 
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which can be written as 
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Using the general relativistic upper bound for the com- 
pactness limit of the star, M/R < 4/9, it is found from 
Eq. (|5U)) that 6(a) < 0. The fact that the parameter C is 
reduced by five-dimensional effects has important conse- 
quences. As it is shown in figure 1, the decrease in C due 
to the CFM exterior solution, unlike the case 8(a) > 
given in Eq. (|3"6")l for the DMPR exterior solution, pro- 
duces a decrese in pressure. Finally, it is easy to see that 
both expressions given by Eq. ([36]) and Eq. (|49| can be 
written in a concise way, as shown below 



8(a) 



1 



3AR 2 (1 + C R 2 )' 2 



y/A(l + Cpfl 2 ) 3 - 1 
y/A(l + C R 2 ) 3 + 1 



f(Co), 

(51) 



where n = for the DMPR exterior solution and n 
for the CFM exterior solution. 



V. CONCLUSIONS 

In the context of the Randall- Sundrum brancworld, 
the consequences of bulk gravity on compact stellar dis- 
tributions was studied through a new brancworld solu- 
tion to Einstein's field equations. This solution was con- 
structed from a well known spherically symmetric stellar 
solution for a perfect fluid in general relativity, namely, 
the Heintzmann solution. In order to generate the 
braneworld version containing the anisotropic effects nec- 
essary for realistic stellar models, the minimal geomet- 
ric deformation approach was used to break the perfect 
fluid solution represented by the Heintzmann solution. 
Hence some important features of the five-dimensional 
consequences on physical variables inside compact dis- 
tributions was clarified. In order to elucidate the spe- 
cific role played by each exterior Weyl function, namely 
the dark radiation U + and dark pressure "P + on compact 
stellar systems, the braneworld solution generated within 
the compact distribution was matched with two differ- 
ent exterior solutions, namely, the Dadhich, Maartcns, 
Papadopoulos and Rezania solution (DMPR) and the 
Casadio, Fabbri and Mazzacurati solution (CFM). Both 
exterior solutions showed different effects on the stellar 
system. It was found that the DMPR solution produces 



an increase in the pressure and that the CFM reduces it. 
The fact that the CFM solution has no dark radiation, 
which is the main difference with the DMPR solution, 
seems to be the most likely cause of the opposite effects 
of these two exterior solutions on the stellar distribution. 
The physically acceptable solution developed in this 
paper represents the point of view of a brane observer, 
hence it is not known whether the bulk eventually con- 
structed will be free of singularities. Despite the above, 
the important thing here is the fact that it was found a 
strong evidence showing that the exterior dark radiation 
IA + always increases both the pressure and the compact- 
ness of the stellar structures, and that the exterior dark 
pressure V + always reduces them. If this is a general fea- 
ture, it would mean that both exterior Weyl functions IA + 
and V + have well defined consequences on stellar struc- 
ture. Therefore an exterior solution with IA + = and 
V + =/= surrounding a stellar distribution might be seen 
as an environment whose physical effects on the stellar 
structure are such that it can be considered as a region 
with negative effective pressure. This is certainly an in- 
teresting case which deserves further investigation. 
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